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ABSTRACT 

We propose observations of the molecular gas distribution during the era of reionization. At z ~ 6 to 8, 
'^CO(y = 1 -0) line intensity results in a mean brightness temperature of about 0.5 fiK with rms fluctuations of 
0.1 i-iK at 1 to 10 Mpc spatial scales, corresponding to 30 arcminute angular scales. This intensity fluctuations 
can be mapped with an interferometer, similar to existing and planned 21 -cm background experiments, but 
operating at ^ 12 to 17 GHz. We discuss the feasibility of detecting the cross-correlation between HI and 
CO molecular gas since such a cross-correlation has the advantage that it will be independent of systematics 
and most foregrounds in each of the 21 -cm and CO(l -0) line experiments. Additional instruments tuned to 
higher-order transitions of the CO molecule or an instrument operating with high spectral resolution at mm- 
wavelengths targeting 158 /im CII could further improve the reionization studies with molecular gas. The 
combined 21 -cm and CO line observations has the potential to establish the relative distribution of gas in the 
inter-galactic medium and molecular gas that are clumped in individual first-light galaxies that are closely 
connected to the formation of massive stars in these galaxies. 



Subject headings: cosmology: theory — diffuse radiation 



1. INTRODUCTION 

Recently there has been a great deal of interest in the evolu- 
tion of the intergalactic medium at high redshifts (z > 6) and 
specifically, the reionization history of the Universe. Most of 
this interest is sparked by existing multi-wavelength observa- 
tions that suggest a complex reionization history, with a strong 
possibility for extended and highly inhomogeneous period of 
reionization. Observations of the 21 -cm spin-flip line of neu- 
tral hydrogen are currently considered to be one of the most 
promising probes of the epoch of reionization (e.g., Madau et 
al. 1997; Loeb & Zaldarriaga 2004; Gnedin & Shaver 2004). 
Given the line emission, leading to a frequency selection for 
observations, the 21 cm data provides a tomographic view of 
the reionization (Santos et al. 2005; Furlanetto et al. 2004a). 
It is also a useful cosmological probe (Mao et al. 2008; Santos 
& Cooray 2006; McQuinn et al. 2006; Bowman et al. 2007). 

Here we propose spectral line intensity mapping associated 
with molecular gas during the era of reionization as an ad- 
ditional probe during the formation of first-light galaxies. In 
particular, we study the 7 = 1 to transition of the '^CO (car- 
bon monoxide), which has a rest- wavelength of 2.61 mm. Our 
suggestion extends the initial work on this topic by Righi et 
al. (2008; see also Basu et al. 2004). Instead of the angular 
power spectrum, which is dominated by CMB anisotropics 
at tens of arcminute scales, we propose spectral line inten- 
sity mapping leading to the three-dimensional power spec- 
trum (see also Visbal & Loeb 2010). An additional probe 
of reionization is useful since the first-generation 21 -cm ex- 
periments are expected to be noise-dominated and the obser- 
vations could, in principle, be affected by various systematics 
and residual foregrounds. A second probe of the reionization 
epoch providing three-dimensional information can be cross- 
correlated with the 21 -cm spectral line measurements to im- 
prove the overall understanding of reionization. 
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In terms of potential spectral lines of interest, the line inten- 
sity of the / = 1 to transition of '^CO has been used as a way 
to establish the gas mass of star-forming galaxies both at low 
(e.g., Downes & Solomon 1998) and medium (e.g., Greve et 
al. 2005) redshifts. Individual CO line observations now exist 
out to z = 5.3 (e.g., Riechers et al. 2010), which is significant 
given that the highest 21 -cm line emission from a galaxy cur- 
rently measured is at z of 0.24 (Lah et al. 2007; see Chang et 
al. 2010 for 21 -cm intensity variations of the kind we propose 
at z ^ 0.9). In addition to the CO lines, we also find that the 
fine structure CII line, with observations in the mm range of 
the spectrum, could be a strong probe of reionizaion. 

This Letter is organized as follows: in the next Section, we 
outline the calculation related to the strength of the expected 
signal from reionization. Section 3 presents results related to 
a potential CO brightness temperature fluctuation experiment 
that target 1-0 transition, and Section 4 discusses the cross- 
correlation with the 21 -cm background. 

2. THE CO( 1-0) SIGNAL 

For a galaxy at a given redshift the CO(l-O) luminosity is 
taken to be related to the H2 content and the metallicity of the 
gas (Obreschkow et al. 2009b). A key input here is the re- 
lation between CO(l-O) luminosity, Leo, and the halo mass. 
While such a relation is expected it also has a large scatter 
coming from differences in the gas content, at a fixed halo 
mass, due to variations in the halo ages at a fixed mass; The 
gas in older halos is more likely to be depleted more than the 
younger ones with older halos containing more stars. How- 
ever, for sufficiently large volumes, we will be averaging over 
many galaxies so that, so the the mean relation and its scatter 
can be quantified. This is similar to the luminosity-halo mass 
relation for optical light which has allowed statistics such as 
the conditional luminosity functions to be constructed and 
improve the halo modeling beyond the simple halo occupa- 
tion number (e.g., Yang et al. 2004; Cooray & Milosavljevic 
2005b; Cooray 2006). 

Here we make use of the numerical simulations of 
Obreschkow et al. (2009d) that are available as part of the 
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Figure 1. Left: The Lco(^hiilo) relation for CO(l-O) luminosity as a function of the halo mass Mhaio at z = 7. The hnes shows the mean relation (thick center 
line) and itlcr relation (thin lines) given the scatter (see text for details). The dots show the mean of the scatter when binned to 150 logarithmic intervals in halo 
mass. Right: The relation neutral hydrogen mass content in individual galaxies Mm vs. the halo mass Mhaio at z = 7. 

SKA Simulated Skies website'*. This simulation is based on 
the galaxy catalog from De Lucia & Blaizot (2007) which 
post-processed the Millennium dark matter simulation with 
semi-analytical models of galaxy formation. It captures the 
gas astrophysics, especially HI and H2 and allow for the cos- 
mic evolution of these two phases of the cold hydrogen gas 
(see, Obreschkow et al. 2009b; 2009c for details). They also 
constructed a virtual sky field with HI and CO line luminosi- 
ties for each of the galaxies, calculated from a combination 
of the neutral and molecular hydrogen gas masses and the 
metallicity of each galaxy. Using the online tools we made 
a query with a four degree observing cone containing more 
than 5 X 10'' galaxies at each of z = 7. 

The Lco{M) from every halo at z 7 from this simulation 
is shown in the left panel of Figure [T] with the shade of gray 
darker in the regions where the halo density is higher Typ- 
ically we find one galaxy per halo at these redshifts. The 
lines show the result of averaging the total Leo from the 
halos for each of 150 logarithmic bins of halo mass. Mo- 
tivated by parameterization of optical luminosity and halo 
mass (Cooray & Milosavljevic 2005a), we use a relation of 
the form Lco(M) = Lo(M/Mcf(l +M/M,)-'' to describe the 
mean relation and determine the four free parameters A, b, c 
and d. At z = 6,7, and 8, these parameters take the values of 
Lo = 4.3 X 10^6.2 X 10^4.0 x lO*" L©, b = 2.4,2.6,2.8, = 
3.5 X 10",3.0x 10",2.0x 10" M©, and c/ = 2.8, 3.4, 3.3, re- 
spectively. 

Using the above relation between Leo and halo mass M as 
a function of redshift, we can write the mean intensity of the 
CO( 1-0) line as 
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Figure 2. The bias and brightness temperature of CO line for z = 6, 7 and 
8. The 68.3% C.L. (Icr) are also shown in green dashed lines. The bias has 
no additional uncertainty as the uncertainty in Lco{M) relation is a change in 
the overall amphtude at a given mass and not a shape change (see. eq. 2). 



^CO = 



(1) 



where we take M^i„ = lQ^{MQ/h) as the minimum mass of 
the dark matter halos that can host galaxies (to be consistent 
with mass scale of atomic Hydrogen cooling or a virial tem- 
perature of lO'^K (Loeb & Barkana 2001), dn/dM is the halo 
mass function (Sheth & Tormen 1999), Di is the luminosity 
distance, and Dp, is the comoving angular diameter distance. 

http://s-cubed.physics.ox.ac.uk 



In above y{z) = dx/di^ = Aco(l +1}^ /H{z), where x is the co- 
moving distance, v is the observed frequency, Aco = 2.6 mm 
is the rest frame wavelength of the CO(l-O) line. 

The signal is expected to contain spatial variations around 
the average intensity due to poisson fluctuations in the number 
of halos (shot noise) and correlations with the underlying dark 
matter density field which will enhance the fluctuations by 
a factor of (\ + b{z,M)5), with 6 being the density contrast 
of this dark matter field. The intensity of the CO(l-O) line 
can be written as a function of the spatial location /co(x) = 
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Figure 3. Left panel: the dimensionless auto power spectrum and the shot noise power spectrum for CO Hne at j = 7. The 68.3% C.L. (Icr) is shown in green 
thin lines. Middle panel: the dimensionless auto power spectrum and the shot noise power spectrum for 21-cm emission at z = T- The auto and shot noise power 
spectrum of the galaxies and A*°' are so small compared with A[p^, that AjJ' is dominated by Aj^'''. The blue lines with error bars make use of noise 
power spectrum (shown by the blue dashed line) coiresponding to the design of LOFAR. Right panel: the dimensionless cross power spectrum and shot noise 
power spectrum for CO line and 21-cm emission at z = 7. The blue solid line is the error from LOFAR. The 68.3% C.L. (Itr) is shown in red thin lines. The 
long-dashed magenta lines in the middle and right panels are the results from the simulation of 21-cin in Santos et al. (2010) for comparison. 



^co[l +^co^(x)] with bias bco given by 



^co(z) = 



Jm 



dM^Lco 



(2) 



where b(z,M) is the halo bias (Sheth & Tormen 1999). 

In Figure 2, we show both the mean signal and the halo 
bias as a function of the redshift during reionization. For sim- 
plicity, we have converted the mean intensity of the signal to 
Rayleigh-Jeans brightness temperature. At z ^ 6, the signal is 
at the level of 0.2 to 1 .5 fiK, with the large range coming from 
the scatter in the Lco{M) relation. 

Since the mean intensity is challenging to measure as it re- 
quires absolute measurements we discuss the feasibility of 
measuring anisotropics of the CO line intensity. The total 
power spectrum of the CO(l-O) line involves two contribu- 
tions, the clustering term and the shot noise. Following the 
mean intensity and bias given above, the clustering power 
spectrum is 

P^^%z,k) = fiobloPss(z,k). (3) 

Here fco is the mean temperature of CO line, Pgs = P\in(z,k) 
where P\in is the linear power spectrum of the dark matter. 
The shot niose power spectrum, due to discretization of the 
galaxies, is 
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3. CROSS POWER SPECTRUM OF CO LINE AND 21-CM EMISSION 

In order to calculate the expected strength of the correla- 
tion between 21-cm and CO line intensities, we briefly discuss 
the brightness temperature variations in the 21-cm signal. We 
consider two signals, the usual signal coming from neutral gas 
in the intergalactic medium that is not directly correlated with 
molecular gas in galaxies, and the remaining neutral gas in 
individual galaxies. During the era of reionization the 21-cm 
brightness temperature fluctuations are expected to be domi- 
nated by the IGM. At lower redshifts. While the IGM signal 
is zero, a small 21-cm contribution still remains due to the 
neutral gas in individual galaxies. 

The difference of the spatially averaged brightness temper- 
ature of 2 1 -cm emission in the IGM and the CMB temperature 



at redshift z (Mao et al. 2008) is 
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Here is the mean neutral hydrogen fraction, that we as- 
sume xh = 1 -Xi where Xi is the mean ionized fraction, and 
we set xu = 03 and 0.7 at z = 7 and 8 respectively (Mao et al. 
2008). The is the averaged spin temperature of the IGM, 
and we assume ^ Tcmb here. The parameter c(z) (defined 
in Santos et al. 2008) is 
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where fib is the baryon density parameter. 

We obtain the power spectrum of the 2 1 -cm fluctuations in 
the IGM using the fitting formula of Mao et al. (2008) at 
several redshifts during reionization. For comparison we also 
consider the reionization simulations of Santos et al. (2010). 
These simulations have a smaller bubble size on average than 
the models of Mao et al. (2008). As we will discuss later, 
while two models we consider give rise to the same 21-cm 
power spectrum, they different significantly in terms of the 
cross-correlation of CO and 21-cm brightness temperature 
fluctuations. 

Similar to the IGM calculation the mean brightness temper- 
ature of 21-cm emission from the galaxies can be calculated 
through 

/ X Ph(z) , 



Xpkiz) 



(7) 



where X = 0.76 is the hydrogen mass fraction, Phiz) = fib(l + 
z)^Pc is the baryon density at z and p^ = 2.7752 x lO" M^/h 
(Mpc//i)-l 

We assume that the total hydrogen mass in a given pixel of 
the experiment only depends on the halo mass so that the hy- 
drogen mass density pH is given by pn = /^___ dMj^MniiM). 
To obtain Mhi(M) we used the same procedure as before, 
querying the simulation from Obreschkow et al. (2009d) and 
adding the Mui from the galaxies in each halo. The result- 
ing scatter plot is shown in the right panel of Figure [T] to- 
gether with the average in each mass bin (blue dots). The 
red Hne shows a fitting function to the average values using 
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Mm = l(yixi°SioW+fi, where M is the halo mass, A = 0.75, 0.74 
and 0.75, and B = 0.17, 0.23 and 0.07 for z = 6, 7 and 8 re- 
spectively. 

Following the description related to the CO line intensity, 
the mean 21 -cm brightness temperature is taken to have spa- 
tial variations IJp(x) = T^(l+bii6(x)), where is the bias of 
galaxies containing neutral hydrogen 
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The clustering power spectrum of the 21 -cm emission from 
galaxies is then P^(z,k) = (T^fbl^Pss(z,k). Again, the shot 
noise power spectrum for 21 -cm emission is 

2 

(9) 
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The total power spectrum of 21 -cm emission follows from a 
sum of shot-noise and clustering terms (see middle panel in 
figure O. 

The cross power spectrum between 21 -cm and CO line in- 
tensities is jj = ^'ccfn+'Pco H+-fco H' where we again sep- 
arate the contribution from neutral hydrogen in the IGM and 
neutral hydrogen remaining in individual galaxies. Based on 
the discussion related to 21 -cm fluctuations, the cross power 
spectrum of the CO line and 21 -cm emission in the IGM is 

Pmfliz.k) = fcobcocizm -x,)Pss-x,Psjl (10) 

where Ps,s(z,k) is the cross power spectrum for the ionized 
fraction and the dark matter The IGM gas is correlated with 
CO lines as they both trace the same underlying dark matter 
density field. 

The cross power spectrum of the CO line and 21 -cm emis- 
sion in individual galaxies is 



Pco,}i(z,k) = TcoT°bcobHPssiz,k), 



(11) 



while the shot noise power spectrum of the CO Une and 2 1 -cm 

emission in the galaxies is 
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4. DISCUSSION 



In Figure 3, left panel, we show the brightness tempera- 
ture power spectrum for the CO line at z = 7, while in the 
middle panel we show the 21 -cm brightness temperature fluc- 
tuations at the same redshift and in the right panel, we show 
the cross-correlation between the two. In the case of 21 -cm 
fluctuations, the solid lines show the prediction based on the 
models of Mao et al. (2008), while in dashed lines we show 
the prediction from the semi-numerical simulation of Santos 
et al. (2010). 

To discuss the possibility of measuring the CO signal and 
the CO-21 cm cross-correlation we make use of a calcula- 
tion similar to the one used in 21 -cm interferometers. Fol- 
lowing Santos et al. (2010), the noise power spectrum, P^ 
for a given experimental "pixel" in Fourier space is given 
by PN(k,e) = Dly(\%%)/[A%n (DA^sin(0)/2^)] where A, 
is the collecting area of one element of the interferometer 
(which could be a station), fo is the total observation time 
and the function n() captures the baseline density distribu- 
tion on the plane perpendicular to the line of sight, assum- 
ing already it is rotationally invariant (k is the moduli of the 



wave mode k and 9 is the angle between k and the line of 
sight). Typically this baseline distribution will be more con- 
centrated around the core but for interferometers with high 
filling factors we can make the simplification that the distri- 
bution is constant so that the required function for above is 
«() = X^N^ /nOl^g^ with Na as the number of elements of the 
interferometer and D,„ax the maximum baseline of the distri- 
bution. Note that the error in the measured power spectrum 
at a given k is just (Ps + Pn) /y/{N,n), where Ps is the signal 
power spectrum and A',,, is the number of modes contributing 
to a given binned measurement. In order to count these num- 
ber of modes, we consider a grid in k space with a resolution 
given by (27r)-'/(yr^ -FoV • B), which is set by the volume of 
the experiment (B is the bandwidth used in the analysis and 
FoV is the field of view). 

Typically we will be interested on linear scales (k < 1.0 
h/Mpc) corresponding to a few arcminutes. Given the fre- 
quencies involved, the need for large FoV, and high sensitivity 
to /iK brightness temperature fluctuations leads to an exper- 
imental setup that involves a total collecting area Atot = 385 
m^, bandwidth B = I GHz with spectral resolution Siy = 30 
MHz and interferometer spacings between 0.7 and 25 m. The 
number of elements is 1000 with each having a receiver with 
system temperature of 20K. In the left panel of Figure 3, we 
assume a total integration time of 3000 hours when calculat- 
ing the noise power spectrum. For 21 -cm noise calculations 
we make use of a setup similar to LOFAR, though we note that 
our calculations are not too different for an experiment like 
MWA. For experimental parameters related to LOFAR, we re- 
fer the reader to Harker et al. (2010). The values (at 150 MHz) 
are Atot = 5 x 10"* m^, T^y^ = 490K (instrument noise tempera- 
ture), FOV=25 deg2, bandwidth of 8 MHz with Si^ = 0.5 MHz 
and Dmax = 2000 m. We take Tint = 1000 hours. 

For the proposed experimental setup the brightness temper- 
ature fluctuations can be detected with a cumulative signal- 
to-noise ratio that leads to more than 40cr confidence. How- 
ever, CO(l-O) mapping could be impacted by foregrounds, 
including spectral lines associated with other molecules; In 
the case of other transitions of the CO lines, the contamina- 
tion for CO(l-O) measurements at z ~ 6 to 8 is likely to be 
minimal; the CO(2-l) line for example must originate from 
z ^ 13 to 17. The most significant contaminations will likely 
come from molecular transitions that have rest wavelengths 
long ward of the CO(l-O) line and thus coming from lower 
redshifts. To avoid such contaminations and to improve the 
overall study of reionization we suggest a cross-correlation 
between 21 -cm brightness temperature fluctuations and the 
CO intensities. 

For the cross-correlation, we assume that observations will 
be conducted in overlapping areas on the sky by both the CO 
and 21 -cm experiment. This is a necessary feature of cross- 
correlation study leading to some coordination between two 
experiments. We reduce both datasets to a common, and a 
lower resolution cube in frequency and count the total num- 
ber of common modes contributing to a given bin in k space 
using the same method for noise calculation as before for each 
of the experiments. Note that the error in a given pixel in k 

space can be written as ^ (P^q jj + Pco^'h) /N,,,, where N,,, is 

the number of modes. In Figure 3 right label we show the 
cross-correlation power spectrum at z = 7 and the expected 
binned errors by using LOFAR and the CO experiment with 
parameters outlined in above. The overall uncertainties are 
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dominated by the 21 -cm observations. For reference, we also 
show the prediction of the model by Mao et al. (2008) and 
from the numerical simulation of Santos et al. (2010). While 
the two give a similar 2 1 -cm fluctuation power spectrum, they 
differ significantly in terms of the cross-correlation with the 
CO emission. The difference is primarily related to the av- 
erage scale of the reionized bubble sizes and we find that the 
cross-correlation such as the one we propose here is more sen- 
sitive to astrophysics during reionization than with 21 -cm data 
alone. 

While our calculation and results are related to the CO(l- 
0) transition, our method can be extended to other transitions 
and other molecules as well as long as a correction is made 
for the ratio of line intensity between higher J transitions of 
CO and the 1-0 line. Existing observations, at low to moderate 
redshifts, suggest a value of about 0.6 for 7 = 2 to 1 luminosity 
when compared to 1-0 luminosity. The 2-1 observations will 
be over the frequency range of 25.5 to 32 GHz for z ^ 6 to 8. 

The fine structure line Cll has been detected in several 
of the high-z galaxies and existing studies indicate that 
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i'Cii/i'CO ^ 4100. With a rest wavelength at 158 /im, obser- 
vations are required over the frequency range of 210 to 270 
GHz (1.1 to 1.4mm). The rms fluctuations are expected to 
be at the level of 10^ to 10^ Jy/sr at z 7 to 8 at tens of 
arcminute angular scales. To probe reionization fluctuations 
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areas of order a few tens square degrees. In a future paper we 
hope to return to the scientific case of a potential Cll experi- 
ment, since such an observation can probe multiple redshifts 
with both CII and high- J CO lines. 
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